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Rainfall from ice-free cumulus clouds requires collisions of large numbers of microscopic droplets 
to create every raindrop. The onset of rain showers can be surprisingly rapid, much faster than the 
mean time required for a single collision. Large-deviation theory is used to explain this observation. 
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The dynamics of the onset of rainfall from ice-free 
(‘warm’) cumulus clouds is poorly understood [H-Q- 
A rain drop grows by collisions of microscopic water 
droplets. A large number of microscopic droplets must 
combine to make one rain drop: the volume increase is a 
factor of approximately one million. The collision rates 
in the early stages of the growth process are low (typ¬ 
ically of order one collision per hour). Given the large 
number of collisions which must occur, it is very hard 
to understand the observation that rain showers can be 
initiated over relatively short periods, of perhaps twenty 
minutes. 

One possible resolution is a consequence of the large 
number of microscopic droplets which must combine to 
make a raindrop. This implies that only very few drops 
are required to undergo explosive growth, and perhaps 
there are sufficient rare combinations of rapid multiple 
collisions to explain rainfall: this point has previously 
been emphasised by Kostinski and Shaw . A quantita¬ 
tive approach is required to show whether the large num¬ 
ber of collisions required for runaway growth can occur 
with a sufficiently high probability. Because this problem 
involves the analysis of rare events, methods based upon 
large deviation theory [^, are used in this Letter to in¬ 
vestigate the hypothesis that rare combinations of rapid 
collisions trigger showers. It is shown that a rain shower 
can develop over a timescale which is a small fraction of 
the mean timescale for one collision. 

First, consider some observations and estimates P-Q 
which illustrate the difficulties in making a quantitative 
description of rainfall. A convecting cumulus cloud which 
could produce showers may have droplets of mean radius 
radius is oq = 10 ^m, which result from condensation 
onto aerosol nuclei. Raindrops have a much larger size, 
typically 1 mm. The volume of a droplet which becomes 
a raindrop therefore increases by a very large factor, de¬ 
noted by Af, which is typically Af « 10®. The num¬ 
ber density of microscopic droplets is typically of order 
Nq = 2.5 X 10® m“®, which gives a liquid water content, 
expressed as a volume fraction, « AttNoUq/S ~ 10“®. 
The cloud depth may he h = 2 x 10® m and the typical 
vertical velocity of air inside the cloud has magnitude 
U Ri 2ms“^, so that the turnover time for convection 
is approximately Th = 10® s. Rainfall from this type 


of cloud can develop over a timescale of approximately 
20 min k. 10® s. 

Droplets which undergo a geometrical collision (the im¬ 
pact parameter is less than the sum of the radii) might 
not coalesce, because the streamlines of small droplets 
curve around larger ones. In fact, if the Navier-Stokes 
equations were a complete description, droplets would 
never collide, because there would always be a lubricat¬ 
ing film of air between them. The coalescence efficien¬ 
cies e of small droplets are somewhat uncertain, but it 
is widely accepted that they are low for typical cloud 

droplets Bi- If the larger droplet has radius below 

20/rm, it is believed that e < 0.1, and that for radius 
10/rm, £ < 0.03 0. For droplets of size a = 50/rm 
colliding with droplets of size a = 10 /rm, however, the 
efficiencies are expected to be close to unity [l|, . 

Collisions between droplets settling at a different rate 
yield a very small collision rate. The Stokes law for the 
drag on a sphere at low Reynolds number indicates that 
the gravitational settling rate is 

2 2 Pi g , , 

v = Tpg = Ka , K=- - ( 1 ) 
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where Tp is the response time characterising the Stokes 
drag on a droplet, p\ is the density of liquid water, and 
Pg and 1 / are, respectively, the density and kinematic vis¬ 
cosity of air. Inserting values for air and water at 5°C 
gives K Rs 1.4 X 10®m“^s“^, so that when qq = 10pm the 
terminal velocity is w 1.4 x 10“^ ms“^ and the response 
time is 1.4 x 10“® s. The collision rate of a drop of radius 
oi with a gas of droplets of radius oq is 

Ri = neNo{ao + ai)‘^K{al - al) ( 2 ) 

where e is the collision efficiency. Setting oi —oq = 2.5 pm 
and £ = 0.03 in addition to the parameters defined above 
gives i?i Ri 2 X 10“®s“®. The rate of coalescence of typical 
sized water droplets due to collisions is therefore very 
small. 

Cumulus clouds are turbulent because of convective 
instability. Saffman and Turner Q investigated the 
role of turbulence in facilitating collisions between water 
droplets. In the case of very small droplets, the collision 
rate due to turbulence is a consequence of shearing mo¬ 
tion. The shear rate of small-scale motions in turbulence 
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is the inverse of the Kolmogorov timescale, tk = a/J^/ e, 
where e is the rate of dissipation per unit mass. Ac¬ 
cording to the Saffman-Turner model, shear induces a 
collision speed of order oq/tk. They argue that the cor¬ 
responding collision rate is 


.Rturb — 


r&K Noe{2a)^ 
Tk 


( 3 ) 


For the parameters of the cloud model, the rate of dissi- 

’ m^s“^, giving tk ~ 70 ms, 


pation is e ^ C/^/t/j ss 2x 10 3 


which gives i?turb ~ 10“° s“^, which is negligible. The 
effects of turbulence are dramatically increased when the 
effects of droplet inertia are significant: this was noticed 
in numerical experiments by Sundaram and Collins Q, 
who ascribed the effect to a clustering effect termed ‘pref¬ 
erential concentration’ Q. More recent work has pro¬ 
posed an alternative mechanism, which has been termed 
the ‘sling effect’ [l^, and which has been explained in 
terms of the existence of caustics in the velocity field of 
the droplets Uj. Inertial effects are measured by the 
Stokes number, St = Tp/rx. Recent numerical studies 
[l^ (see also 13 1 ) show that the collision rate is greatly 
enhanced by effects due to caustics for St > 0.3, equa¬ 
tion dHI) is a good estimate when St ^ 1. While it is 
in principle possible for turbulence to be responsible for 
an enhanced collision rate of water droplets due to iner¬ 
tial effects, the parameters of the cloud model discussed 
above yield St Ri 2 x 10“^, where there is no signihcant 
enhancement. While there is a consensus that turbulence 
is important for the formation of rain showers 0 , tur¬ 
bulent enhancement of collision rates does not appear to 
be sufficient. 

Now consider how to model the onset of a shower. 
It has already been remarked that showers occur on 
a timescale which may be smaller than the typical 
timescale for one collision. It is, therefore, reasonable 
to assume that the runaway droplets are falling through 
a background of droplets which have not yet coalesced. 
As a runaway droplet falls it collides with a large num¬ 
ber M of small droplets of size oq. The time between 
successive collisions may be assumed to be independent 
Poisson processes. If the time between the collision with 
index n and the previous collision is the time for a 
droplet to experience runaway growth is 


At 


( 4 ) 


i=l 


where the are independent random variables with a 
Poisson distribution 


Pn(^n) — 7?n exp( Rn^n) ■ (5) 

The problem is to determine the statistics of T in the 
limit as A/" —>• oo. The rates for successive collisions in¬ 
crease as the size of the falling drop grows. Because all of 


the collision rates scale in the same way as a function 
of the droplet size oq and the number density Nq, write 

= i?i/(n) . ( 6 ) 

Here Ri depends upon the properties of the cloud but 
the function f{n) is the same for all clouds. In order to 
identify the form of /(n), consider the rate of collision of 
a large droplet resulting from n previous collisions with 
a gas of small droplets of radius oq. The radius of the 
large droplet is an = n^^^ao- When n is large it may 
assumed that the collision efficiency is e ss I and an 2 > og 
so that Rn irNonaf^ oc which suggests setting 

/(n) = However during the early stages of droplet 

growth, the collision efficiency for the first few collisions 
is small, but increases rapidly with n. In what follows 
f{n) is assumed to be a power-law 

Rn = Rm"' ■ (7) 


If the collision efficiency of droplets were unity, it would 
be appropriate to set 7 = 4/3. Because the collision ef¬ 
ficiency of droplets at the crucial initial stage of their 
growth is small, the collision rate increases more rapidly 
as the size of the falling droplet increases. When the 
droplets are between 10 fim and 50 /rm it is reasonable 
to model the product of the collision rate and the colli¬ 
sion efficiency as being proportional to a®, that is to n^, 
where n is the number of collisions [^. In other cases, 
such as solid precipitation (snow), other values of 7 may 
be appropriate. In the following 7 is left as an adjustable 
parameter, but special consideration is given to 7 = 2 , 
because it gives a good approximation to terrestrial rain¬ 
fall, and to 7 = 4/3, because this may be a good ap¬ 
proximation for atmospheres on other planets where the 
collision efficiency might not limit the rate of coalescence. 

It is necessary to determine the probability density for 
the time T being a very small fraction of its mean value, 
(T). Inspired by large deviation theory @ , the proba¬ 
bility density of T is written in an exponential form: 

PiT) = exp[-J(T)] , ^ ■ (3) 

When f{n) is given by ([7]), the mean time for explosive 
growth converges as Af ^ 00 when 7 > 1 : 




( 9 ) 


where ^ is the Riemann zeta function. The function J (t) 
in dH) is often termed the entropy in texts on large devia¬ 
tion theory. It will be necessary to determine the entropy 
function J(t) from the rate function f{n). 

After a drop has grown to a size where it is much 
larger than the typical droplets, and where the colli¬ 
sion efficiency is approximately unity, it falls rapidly and 
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collects other droplets in its path. Consider a drop of 
size oi falling through a ‘gas’ of much smaller droplets, 
with liquid volume fraction $i. The larger drop falls 
with velocity v = Kof and grows in volume at a rate 
TTea\^\v = 47 rafdi, where di is the rate of increase of 
the drop radius. The rate of increase of the radius of the 
‘collector’ drop as a function of the distance x through 
which it has fallen is 


da e$i 
dx 4 


( 10 ) 


Note that this expression is valid whether or not the ter¬ 
minal velocity is given by the small Reynolds number 
approximation, ©■ In the case of droplets which reach 
a radius of approximately 1 mm, the collision efficiency e 
is close to unity throughout most of the fall. The droplet 
radius after falling through a cloud of depth h is therefore 
a{h) ~ $i/i/4. It will be assumed that the most relevant 
collector drops are those that started at the top of the 
cloud, so that the volumetric growth factor is 


Using the representative values given above gives J\f ss 
10 ^ 

The rate of change of the liquid water content of a 
cloud due to the runaway growth of droplets is 


^ . ( 12 ) 

Note that the growth factor J\f and the probability den¬ 
sity for runaway growth after time t are both functions 
of 4)i, but if the objective is to understand the onset of 
a rain shower it suffices to evaluate these quantities with 
the initial value $ 1 ( 0 ). Using ([5]) for P{t), the condition 
for the timescale t* where there is a significant reduction 
in <i>i(t) is Af ejq>[—J{t*/(T))] = 1, that is 


t*=T*{T), J(t*)= In A/". (13) 


The droplet volume growth factor Af was estimated in 
equation m- To determine determine the solution of 
(fT3l) for t *, it is necessary to determine the entropy func¬ 
tion J (t) for the random sum defined by dH) and dU . 

To compute J(T) consider a cumulent generating func¬ 
tion A(fe), defined by writing 

nOO 

exp[—A(A:)] = (exp(—/cT)) = / dT P(T) exp(—AT) . 

Jo 

(14) 

Because the tn are independent, with a distribution given 
by dS]): 


A(fc) = ln(exp(-fct„)) = XI ^ ) 

n=l n=l A nj 


(15) 


Now consider how to obtain P{T) from X{k). Noting that 
exp[—A(fc)] is the Laplace transform of P{T), application 
of the Bromwich integral formula for inversion gives 

^ pC-\-ioo 

P{T) = dz exp[zT - A(z)] (16) 

27ri Jc-ioo 


where C > —Ri. The integral is dominated by contri¬ 
butions from the neighbourhood of a saddle at z = k*, 
where 


M 

^ = E 


Rn 


(17) 


which is to be solved for k* given a value of T. The 
probability density P(T) is then approximated by 


where t = T/ (T) and S{k) is the magnitude of the second 
derivative of the exponent in (1161) . Equations (HD and 
dH cannot be solved exactly and explicitly for J(t). 
Consider how to write down a parametric representation 
of J(t) using a scaled variable, k, defined by k = k*/Ri. 
The dimensionless time for raindrop formation is 


AT 


"(«) = X 
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n—1 


K + n'^ 


■ AT 

x^ 




(19) 


and the entropy function is 


AT 


AT 


'^(^) = X (1 + - X 




n—1 


n—1 


1 + Kn 


( 20 ) 


Figured!] shows the distribution of r = T/{T) for the case 
7 = 4/3, with Af = 4000 and i?i = 1, comparing the re¬ 
sults of simulation of (jd]) , the Bromwich integral (I16|) , the 
saddle-point approximation, equations dUD, dUD, (EOD, 
which are all in excellent agreement. Figure [2] makes a 
similar comparison for the case 7 = 2 , which is most rel¬ 
evant to rain showers. In both cases the entropy function 
increases very rapidly as r —>■ 0 , indicating that the value 
of T* = t* /{T) is quite insensitive to the value of In A/”. 
It is clear from figure [2| that the solution of equation (fT^ 
gives small values of r* when Af is large. Numerical eval¬ 
uation of (|16D with 7 = 2 gives r* « 0.077 when A/” = 10® 
and T* « 0.068 when Af = 10®. Alternatively, in terms 
of (U) = when 7 = 2 , the predicted time for 

onset of a shower is a small fraction of the mean time 
for the first collision: t* Ri 0.128(ti) when Af = 10® and 
t* Ri 0 . 112 (ti) when Af = 10®. 

It is clear from figures [T] and [2| that the entropy in¬ 
creases very rapidly as r —> 0 , and it is desirable to find 
asymptotic behaviour of J{t) when r <C 1. This limit 
corresponds to k ^ 1. In the limit where Af ^ 00 and 
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FIG. 1. Plot of ln[P(r)], for M = 4000 and 7 = 4/3. The 
results of simulation, evaluation of the Bromwich integral, and 
the saddle point approximation are in excellent agreement. 
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FIG. 2. Plot of ln[P(r)], for M = 4000 and 7 = 2. 


At ^ I in equations (HU) and (EHl), the sum in the numer¬ 
ator of (HH) is approximated by an integral: 


r(At) 


C(7) 


dn 


I 


K + •nH 


ip(7) (21) 


with 


F(7) = 


7C(7) 


dx 


xi ^ 


\ + X 


( 22 ) 


A similar approach applied to (EUl) gives 


J(k) = k^C?( 7 ) , G(-f) = -f dx ^’^^\^ ^^ -C( 7 )-F( 7 ) ■ 

7 Jo x~i~^^ 

(23) 

Eliminating At from (I22II an (l23ll shows that in the limit as 
T —A 0 the entropy function has a power-law divergence: 
to leading order 


J(t) ~ G'( 7 )[F( 7 )] T-i r =K{'^)t . (24) 


This indicates that the probability density has a non- 
algebraic singularity as r —A 0: P(t) ^ exp(—AT/r) 
when 7 = 2, or P(r) ^ exp(—AT/r^) when 7 = |. 
Numerical integration gives A'(|) « 0.8799, whereas 
K(2) = 'K^j4((2) = I exactly. 

The conclusion is that rain showers can commence in a 
timescale which is short compared to the mean time for 
the first collision between droplets, with the timescale for 
onset being approximately one eighth of the mean time 
for first collision in the case of the more realistic model 
(7 = 2). Thus large deviation theory has resolved an 
apparent paradox of meteorology, that rain showers can 
start very quickly, on timescale which are short compared 
to typical mean collision times. 

This calculation does not resolve all of the uncertain¬ 
ties about initiation of rain showers. Clouds can exist for 
a long period without producing a rain shower, before 
depositing a large fraction of their water content over a 
short time. Shower activity is associated with convective 
motion in clouds, and it has been suggested that tur¬ 
bulence facilitates collisions. For typical levels of turbu¬ 
lence, however, turbulent enhancement of collisions does 
not appears to be sufficient. It seems as if non-collisional 
mechanisms involving convection must may a role in ini¬ 
tiating the cascade |15| . 
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